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SUMMARY 


Leadiiig  edge  devices  and  slats  are  applicable  to  high  speed  wing 
planforms  with  high  sweepback  or  strakes.  This  includes  aircraft  with 
transonic  manoeuvrability  capability  and  supersonic  cruise  vehicles. 

A  theoretical  approach  within  the  framework  of  conical  slender  wing 
theory  has  been  presented  that  deals  with  thin  delta  wings  with  leading 
edge  slats  including  flow  separations. 

Specimen  calculations  suggest  that  benefits  in  performance  (lift/ 
drag)  are  strongly  dependent  on  the  flow  conditions  prevailing  at  the 
edges  of  the  configuration  (i.e.  whether  the  flow  is  attached  or  separ¬ 
ated)  . 

In  general  lift-drag  benefits  of  10%  -  40%  are  obtainable  with 
suitable  configurations.  These  benefits  are  of  the  same  order  as  those 
indicated  by  experiments. 

Further  work,  both  theoretical  and  experimental  leading  to  practical 
applications  has  been  recommended. 
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INTRODUCTION 


Leading  edge  devices  and  slats  are  applicable  to  a  number  of  high 
speed  wing  planforms  with  high  sweepback  or  strakes  (Fig.  1).  This 
includes  aircraft  with  transonic  manoeuvring  capability  and  supersonic 
cruise  vehicles  not  only  for  civil  or  transport  functions  but  also  for 
strategic  reconnaissance  and  advanced  air-defence  concepts  (e.g.  against 
cruise-missile  carriers  and  penetrating  bombers).  In  the  far  future, 
liquid  hydrogen  fuel,  if  it  is  to  be  used,  makes  more  sense  for  an  air¬ 
craft  with  supersonic  cruise  capability.  Leading  edge  devices  may  also 
be  used  on  the  canard  wing  of  forward-sweep  wing  aircraft  designs. 

The  leading  edge  design  on  aircraft  is  subject  to  compromises 
between  the  conflicting  requirements  of  high  speed  lift/drag  efficiency 
and  low  speed,  and  transonic  manoeuvring  requirements  of 
high  lift.  The  usual  way  of  incorporating  leading  edge  devices  is  design 
the  wing  for  cruise  or  high  speed  with  devices  retracted  and  then  to 
deploy  them  during  the  phase  of  flight  when  high  lift  is  called  for, 
i.e.  at  low  speed  and/or  transonic  manoeuvres. 

Simple  leading  edge  devices  in  the  form  of  flaps  or  variable  camber 
may  be  used,  but  they  allow  only  one  degree  of  movement  and,  therefore,  a 
limited  scope  for  optimisation.  Leading  edge  slats  offer  one  extra 
degree;  the  position  ana  angle  can  both  varied  for  optimised  condition. 

For  conventional  aircraft  with  wings  of  low  sweep-back,  there  is  a 
considerable  amount  of  information  available  on  the  subject  of  leading  edge 
design  with  or  without  devices  (Ref.  1).  Calculation  methods  have  been 
devised  for  treating  2-D  geometries  (Figs  2  and  3).  The  methods  are 
strictly  applicable  only  when  3-D  effects  such  as  those  due  to  wing-tips  or 
fuselage  junctions,  are  small. 

For  wings  of  higher  sweep-back,  however,  there  seems  to  be  only  a  small 
amount  of  information  available,  particularly  with  regard  to  leading  edge 
devices.  An  idea  of  the  order  of  gains  from  leading  edge  devices  can  be 
obtained  by  reference  to  the  work  of  Ray  and  Hollingsworth  (Ref.  2)  on  F-4 
Fighter  aircraft  with  leading  edge  sweep-back  51.4°  (Fig.  4).  They  conclude 
that  incorporation  of  devices  resulted  in  a  sizeable  33%  improvement  on  the  buf¬ 
fet  onset,  and  L/L  performance  gain  of  35%  at  Cg  =  0.8, and  improved  lateral 
directional  characteristics  throughout  the  test  Mach  number  range  of  0.6  to 
0.94.  The  improvements  were  verified  in  a  subsequent  flight  evaluation. 

Goodmanson  and  Gratzer  (Ref.  3)  show  in  Fig. 5,  24%  improvement  in  L/D 
at  Cg  =  0.4,  using  droop  and  slats  on  a  highly  swept -back  wing  (L.E.  sweep- 
back  70°  approximately).  Experience  of  NASA  Langley  on  supercruiser  and 
arrow  shaped  wings  also  suggests  similar  improvements. 

An  idea  of  complexity  of  such  flowfields  can  be  gained  from  Fig.  6 
which  shows  multiple  vortex  separations  that  arise  due  to  the  use  of 
leading  edge  slats.  The  vortex  separations  are,  of  course,  sources  for  non¬ 
linear  lift  and  provided  this  lift  is  distributed  on  forward  facing  surfaces 
we  have  then  a  mechanism  for  reducing  drag. 


J 

In  a  practical  flow  situation,  compressibility,  vortex  wake  and  j 

boundary  layers  on  surfaces  introduce  further  complexity  to  the  flow 
field . 

(4) 

Recently  Cuming  and  Dickison  nave  investigated  flow  around  delta 
wings  with  conical  type  leading  edge  slats.  Their  studies  have  included 
flow  visualisation  using  water  tank,  smoke  and  surface  oil  flow  patterns. 

They  have  also  conducted  force  tests  on  a  selection  of  wing-slat 
geometries.  They  indicate  that  multiple  flow  separations  and  inter¬ 
actions  are  present  and  benefits  in  L/D  of  10%  -  15%  are  obtainable  for 
certain  conical  type  geometries. 
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II. 

II. 1 


PHYSICAL  FLOW  FIELDS  AND  THEORETICAL  MODELS 


QUALITATIVE  FEATURES  OF  FLOW 

It  is  instructive  to  look  at:  the  qualitative  features  of  flow  that 
lead  to  the  possible  use  of  leading  edge  devices,  e.g.  plain  flaps, 
smooth  variable  camber  or  slotted  flaps  (slats).  It  is  convenient  to 
confine  our  attention  to  conical  flows  in  small,  moderate  and  higher 
incidence  ranges.  Conical  streamline  patterns  on  basis  of  Ref. 9  have  been  sketcl 

Small  Incidence  Range 


The  simplest  situation  of  attached  flow  over  the  wing  with  the 
leading  edge  flap  deflected  to  prevent  any  separation  at  a  particular 
incidence  is  illustrated  in  Fig.  7(a).  This  type  of  flow  may  be  calculated 
using,  for  example,  techniques  such  as  conformal  transformation  or  panel 
methods.  The  flow  is  susceptible  to  off-design  variation  and  a  shoulder 
separation  as  shown  in  Fig.  7(b)  may  appear.  This  shoulder  separation  may 
be  avoided  by  using  slots  near  the  leading  edge  as  shown  in  Figs  8(a)  and 
(b).  In  Fig.  8(a)  the  'rolled-up'  wake  of  the  flap  is  shown  near  its 
trailing  edge  whilst  in  Fig.  8(b)  t his  wake  is  allowed  to  extend  further 
over  the  wing.  The  qualitative  features  of  the  Figs  8(a)  and  (b),  particu¬ 
larly  the  latter,  resemble  those  of  Fig.  7.  It  is  difficult  to  lay  down 
any  rules  about  lift/drag  (L/D)  performance  of  the  flows  at  this  stage. 
Intuitively,  we  might  say  that  the  performance  of  the  wing  with  the  slotted 
flap  would  be  worse  than  the  calculated  performance  of  Fig.  7(a),  but 
perhaps  better  than  the  actual  experimental  performance  of  the  shoulder  flap 
(Fig.  7(b)). 

Moderate  Incidence  Range 

At  moderate  incidences  the  flow  separates  at  the  leading  edge  in  spite 
of  flap  deflection  and  the  flowficld  as  illustrated  in  Fig.  9  may  appear. 
The  separations  near  the  shoulder  and  the  leading  edge  may  be  controlled  by 
using  slots,  as  in  Figs  8(a)  and  (1>). 

High  Incidence  Range 

At  higher  incidence  the  primary  vortex  moves  inboard  of  the  shoulder. 

A  simple  model  is  shown  in  Fig.  iO(a)  and  a  more  realistic  one  in  Fig.  10(b) 
with  a  secondary  separation  fixed  near  the  shoulder.  This  secondary  separ¬ 
ation  could  be  removed  by  a  slot  as  in  Figs  11(a)  an  I  (b).  These  slots  are 
in  the  opposite  sense  to  those  of  Figs  8(a)  and  (b).  They  enable  the 
idealized  flow  of  Fig.  10(a)  to  be  realized,  however. 

It  is  evident  from  the  foregoing  that  there  are  a  number  of  areas  in 
which  investigations  for  understanding  the  flow  are  needed,  both  theoret¬ 
ically  and  experimentally . 

An  idea  of  configurations  suitable  for  modelling  by  theory  can  be 
gained  from  Figure  12, 


II. 2  APPLICABLE  DESIGN  METHODS 


As  far  as  the  design  methods  are  concerned,  the  current  state  of  art 
may  be  inferred  from  a  paper  by  Tinoco  and  Yoshihara^)  in  which  they  have 
considered  separated  flow  past  slender  thin  wings  with  leading  edge  droop 
and  vortex  control  tabs.  The  capability  of  the  computer  code  using  panel 
methods  is  limited  to  one  vortex  separation  system  only.  Any  extension  of 
the  panel  method  to  simulate  the  wing  and  leading  edge  devices  together 
will  require  a  very  large  number  of  panels  to  permit  sufficient  resolution; 
to  this  must  be  added  the  panels  required  for  simulation  of  multiple  flow 
separations  from  various  edges  of  the  configuration. 

The  method  of  leading  edge  suction  analogy  is  more  suited  to  estima¬ 
tion  of  overall  planform  effects  and  not  to  detailed  description  of  flow 
that  is  required  for  understanding  of  leading  edge  devices. 

It  is  evident  from  the  foregoing  that  we  have  some  way  to  go  before 
we  can  model  the  flows  of  Fig.  6  in  a  fully  three-dimensional  way.  In  a 
number  of  applications  the  wing  trailing  edge  effect  can  be  small,  e.g.  on 
slender  wings  the  first  two-thirds  of  the  flow  is  nearly  conical.  On  wings 
with  strakes,  the  strake  portion  is  mainly  subject  to  slender  wing  principles. 
Slender  wing  theory  forms  a  necessary  platform  for  any  three-dimensional 
calculations  such  as  by  step-by-step  method.  A  more  exact  method  is  by 
inclusion  of  chordwise  terms  thus  simulating  the  3-D  trailing  edge  effect . 


The  slender-wing  approach  is  versatile  and  it  can  utilize  either  the 
conformal  transformation  method  or  Fourier  type  series  formulation  in  the 


cross-flow  plane  -  or  a  combination  of  both  of  these.  The  theoretical 
approaches  of  Brown  and  Michael^®),  Nangia^7’8^,  Smith^9\  Pullin^10*11^ 
Jonesd-2}  an  form  essential  background  to  the  method  mentioned  in 


this  report. 
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III  THEORETICAL  MODEL 

III.l  Model  Geometry 


We  take  a  set  of  Cartesian  co-ordinates  (x,  y,  z)  with  the  origin  at 
the  apex  of  the  wing  (component  A  of  the  combination)  as  shown  in  Figure  13. 
The  x-axis  coincides  with  the  centre  line  of  the  wing,  the  y-axis  lies  in 
the  plane  of  the  wing  measured  in  the  starboard  direction  and  the  z-axis  is 
perpendicular  to  the  wing.  The  angle  of  attack  measured  at  wing  centre¬ 
line  with  respect  to  the  freestream  V  is  a.  The  leading  edges  of  the  wing 
(component  A  of  the  combination)  are  defined  by  y  =  ±  k  x.  The  perturba¬ 
tion  velocities  are  denoted  by  u,  v,  w  in  the  x,  y,  z  system. 

The  starboard  slat  (component  B  of  the  combination)  is  located  near  the 
starboard  leading  edge  of  wing  A.  The  centre  line  of  the  slat  lies  at 
y  =  cQ.x,  z  =  hQ.x.  The  semi -span  of  the  slat  is  cj  x  and  its  inclination 
to  y-axis  is  $  (measured  positive  upwards).  For  convenience,  we  introduce 
axes  yfi  and  zg  to  measure  dimensions  along  and  normal  to  the  slat  as  shown 
in  Figure  14.  The  leading  edge  of  the  slat  lies  at  y  =  (c  +  cj  cos  <|>)  x 
and  z  =  (hQ  +  ci  sin  $)  x.  The  trailing  edge  of  the  slat  Jies  at 
y  =  (cQ  -  cj  cos  <p)  x  and  z  =  (hQ  -  cj  sin  <J>)  x.  The  perturbation  velocities 
are  denoted  by  Vg  and  wg  along  yg  and  Zg  directions. 

The  port  slat  (component  C  of  the  combination)  is  a  mirror  image  of  the 
starboard  slat  in  the  xz-plane.  We  introduce  axes  yc  and  zc  for  this  slat 
as  shown  in  Figure  14. 

III. 2  Boundary  Conditions 

Boundary  conditions  on  the  surfaces,  edges  and  in  the  field  are  as 
follows: - 

111. 2.1  Surfaces 

The  problem  is  restricted  to  slender-body  theory  and  a  solution  is 
desired  which  satisfies  the  condition  of  zero  normal  velocity  in  each  cross- 
flow  plane  on  the  wing-slat  combination  surfaces 

111. 2. 2  Edges 

At  the  edges  of  the  configuration,  there  is  a  choice  of  conditions 
available  which  depends  on  the  complexity  of  the  model  adopted  as  shown  in 
the  following  table. 
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In  the  present  work,  edge  vortex  sheets  are  .,ut.  er^  loyed  and  the 
choice  is  restricted  to  four  possible  types  of  wing-slat  configurations 
as  shown  in  Figure  15. 

1.  Configuration  A-S-A  Attached  flow  at  wing  and  slat  leading  edges, 
vortex-cut  arran^pgjj^at  slat  trailing  edge. 

2.  Configuration  S-S-A  Attached  flow  at  wing  leading  edge,  vortex-cut 
arrangements  at  the  slat  leading  and  trailing  edges. 

3.  Configuration  A-S-S  Attached  flow  at  wing  leading  edge,  vortex-cut 
arrangements  at  the  slat  leading  and  trailing  edges. 

4.  Configuration  S-S-S  Vortex-cut  arrangements  at  all  edges. 

In  addition,  the  configuration  of  slats  only  is  of  interest  as  a 
reference.  Two  possible  configurations  are  (Figure  16): 

1.  Configuration  S-A  Attached  flow  at  slat  leading  edge,  vortex-cut 
arrangement  at  slat  trailing  edge. 

2.  Configuration  S-S  Vortex-cut  arrangements  at  both  edges  of  the  slat. 

It  is  thought  that  the  more  complex  edge  vortex  sheet  models  can  be 
considered  once  an  insight  has  been  gained  from  the  calculation  of  the 
foregoing  configurations. 

In  the  cross-flow  plane,  the  behaviour  of  trailing  vorticity  near 
the  edges  can  be  summarised  as  follows: 

Wing  (0A  =  cos  1  (-  y/kx>) 


Flow  at  Wing  Leading  Edge 

Leading  Edge  Vortex 

Form  of  Vorticity  at  Edge 

Attached 

0 

x) 

«  '  or  cot  BA 

- T - 2— 

_ 

Separated 

Tl 

1 

■ 

*  k or  sin  26a 

8 


^8 


Starboard  Slat  (8g  =  sin-1  ( ^P/cjX  )  ) 


Flow  at  Edge 

Edge  Vortex 

Form  of  Vorticity  at  the  Edge 

Attached  at  Slat  Leading  Edge 

Separated  at  Slat  Trailing  Edge 

0 

r2 

1 

[ci +  yB  .  0b 

-  or  tan  — 

Jc!  -  yB  2 

Separated  Flow  at  Slat  Leading  Edge 

Separated  Flow  at  Slat  Trailing  Edge 

r3 

r2 

■J 

1  -(  \  or  sin  Q 

\cj  yj  B 

Port  Slat 


The  behaviour  follows  from  symmetry  considerations. 

Ill  2.3  Field 

(a)  The  flow  is  unperturbed  at  infinity. 

(b)  In  the  case  of  flow  separations  idealised  as  vortex  sheets,  it  is  necessary 
to  satisfy  that  they  are  stream  surfaces  of  three-dimensional  flow.  For 
the  vortex-cut  arrangement  we  need  to  satisfy  only  that  the  arrangement  as 
a  whole  carries  no  load. 


(c) 


To  prevent  vorticity  being  carried  through  the  apex  from  one  element  of  the 
wing-slat  combination  to  another,  it  is  required  that  the  total  circulation 
about  each  element  of  the  combination  is  zero.  The  wing  by  virtue  of  its 
symmetry  satisfies  this  condition  automatically.  The  condition  needs  to 
be  imposed  on  slats  only. 


In  this  report,  as  mentioned  earlier,  configurations  with  and  without  flow 
separation  from  edges  are  of  interest.  In  the  most  general  case,  i.e.  configur¬ 
ation  S-S-S,  it  is  convenient  to  define  various  quantities  as  follows;  we 
need  to  refer  to  starboard  slat  in  the  first  place  (Fig.  17). 


] 


Wine 


I 

-I 


Quantity 

Location 

- 1 

Sense  1 

Vortex  Tj (x) 

yi (x),  zi(x) 

Positive  anticlockwise 

Cut  drl<x> 

between  (kx,  0)  and 

dx 

{yi(x)»  Zitx)) 

(  v  \ 

Along  wing  surface 

Positive  anticlockwise 

Trailing  vorticity  $w(j —  ) 

about  x-axis 

Bound  vorticity 

Along  wing  surface  -  follows 
from  i.e.  3?w.  3iw 

Positive  clockwise 
about  y-axis 

3  y  3  x 

_ A 

Quantity 

Location 

Sense 

Vortex  r2(x) 

(y2(x),  z2(x)} 

Positive  anticlockwise 

Cut  dr2<l0 
dx 

between  { (c0-cicos<ji  )x, 
(hQ-cisin<J>)x}  and 
(y2(x),  z2(x)} 

Vortex  r 3(x) 

y3(x),  z3(x) 

Positive  anticlockwise 

Cut  dr3(lt) 
dx 

between  {(c0+cicos<(>)x, 
(hQ+cisin<J>)x}  and 
(y3(x),  z3(x)} 

Trailing  vorticity 

Along  slat  surface 

Positive  anticlockwise 
about  x-axis 

Bound  vorticity 

Along  slat  surface  -  follows 

x  .  36s 

from  oP,  i.e.  _ _  _ _ 

S  3y  3X 

Positive  clockwise 
about  y-axis 

The  quantities  on  the  port  side  of  the  configuration  are  defined  on 
considerations  of  symmetry.  It  is  convenient,  however,  to  denote  the  port 
slat  trailing  vorticity  by  6p  and  bound  vorticity  by  yp. 

As  in  the  Brown  and  Michael  model,  the  vortex-cut  representation  leads  to 
a  small  finite  load  at  the  edge  under  consideration.  Zero  load  at  the  edge 
can  only  be  obtained  with  vortex  sheet  representation  by  Smith's^9)  model. 
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IV 


DEVELOPMENT  OF  BOUNDARY  CONDITIONS  ’ 

IV  1.  Surfaces  j 

We  need  to  consider  the  cross-flow  plane  in  which  the  trailing  vorticities 
6W»  65  >  6p  and  any  or  all  of  the  vortices  Tj,  T2,  T3  are  present.  | 

A  general  point  (xA,  y A,  0)  on  the  wing 
The  upwash  equation  is 


+  u5ss(-*a  .^a  .°)  + 


,bA,o) 


+  4  0)  4  ,0)+  VowvC<  -  O 


■1  '  "  >  - 

*  1 

V  - 

-«  1 

where 

is 

the 

upwash 

induced 

due 

to  the  vorticity  Swon 

the  wing 

ft 

t» 

11 

If 

"  "  "  Ss  on 

the  slat 

»» 

t» 

»» 

ft 

"  "  "  €pon 

the  slat 

If 

a> 

»» 

i 

M 

»» 

>  St(x) 

♦  » 

"  a  pair  of  vortices 

CrCx) 

Geometry  considerations  show  that 

^SrCx*,bA,o)  ,0). 

The  upwash  equation  may  now  be  put  in  more  convenient  form  with  slat 
induced  velocities  measured  in  Vg,  <*Pg  system. 

^W(*A  ,^A  .  °)  +  {  ^  >  *Ba)  +  '^BA  .  ^  ? 

♦  [  %5(7A  ,  h*  ,*B*)  '  ^Bgs(  **.**  ,  '»*<? 

4  +  ^(xA,bA)o)  + 

where  (xA,  yBA,  zBA)  is  the  point  (xA,  yA,  0)  in  the  x,  yjj,  zB  axes  system 
(Fig.  18)  and  (xA,  yj^,  zBA)  is  the  point  (xA,  -  yA,  0)  in  x,  yB,  zfi  axes 
system. 


Slat  B  at  a  general  point  (xb,  ys»  0) 


The  upwashes  are  measured  in  direction  and  we  have  for  the  most 
general  case  with  all  Tj,  r2  and  r3  being  present 

+■  o)+  "■grj(.I,,iB,o')+V^«s  =  o 

where  UJgL  is  the  upwash  induced  due  to  the  vorticity  S*j  on  the  wing 


fSr> 


11  11  11  ft  ft  If  |1 


ft  If  ft  tt  tf  ft  ft 


£*  on  the  slat  B 


(Oi> 


"  on  the  slat  C 

(r=l,2,3)  is  the  upwash  induced  due  to  a  pair  of  vortices 

r*(x)  Xi  S&Lx),  ?,-(*) 

is  angle  of  attack  of  the  slat. 

Geometry  considerations  show  that 

^,(*8'  > °)  **  \s  CXft  9  V  /  )  <*»  Zf  +  \s  (*«,  >SC  ,£S ^ 

where  is  the  sidewash  velocity  induced  due  to  vorticity  &s  on  the  slat 

•^8c>?8t  are  de*;ined  i°  fig*  l9* 

The  upwash  equation  may  now  be  put  in  more  convenient  form  with  wing 
induced  and  vortices  induced  velocities  measured  in  tr;ur  system. 


ze  >  he  )  ?*e)  ~  xs  >  ^*8  > 2m) 


+  “’b 

3 


V^.Ss.ol  + 

+  y  *«.)*•*  -  WV •  *.«  } 

"t  V/'Jvit  «  O  .  .  .  .  ( 2 ) 

where  the  point  (xg,  yAfi,  zAg)  is  the  xg,  yg,  zg  in  the  x,  y,  z  system. 


I 
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IV  2.  Field  Boundary  Conditions  of  Zero  Total  Circulation  about  Slat 


j 


) 


(3) 


If  the  flow  is  attached  about  the  slat  leading  edge  then 
/ 

r2  ♦ 

If  the  flow  is  separated  at  the  slat  leading  edge  then 

^  +  J  &  s )  ^i/8  *  ° 

IV  3.  Zero  Force  on  Vortex-cut  Arrangements 

The  condition  of  zero  force  on  vortex-cut  arrangement  follows  from 
Brown  and  Michael  (Ref.  6)  and  may  be  written  in  the  general  form  for  a 
vortex  at  z.  ,^1%)  joined  by  a  cut  to  an  edge  at  ^<».)  Z  (■*)  as  follows 


(4) 


2  ft 

5. 

2- 

X 


V 


ft 

t  — 

x 


‘Sr  h 


3* 


where  and  are  the  velocities  induced  at  the  vortex  position 

due  to  whole  flowfield  (excluding  the  effect  of  vortex  on 

itself). 

We  may  now  write  Equation  (5)  for  each  edge  of  the  combination 


(5) 


Wing  Leading  Edge 

V 

2  li 

<*>. 

X 

V 

Slat  Trailing  Edge 

*  ft  , 

X 

'v 

t 


(6) 


(7) 


It  is  interesting  to  note  that  for  the  particular  case  of  $  =  0  and 
vanishing  a,  we  have 

_.  .  -*0  ,  -*» 

Therefore  '  1 


-  13 


Slat  Leading  Edge 

^  +  Q>  +  q  cwj 

-  V  T 

2  Is  .  •£*  +  K0  +  q  f 

x-  '  v 


-  14  - 


V  SOLUTION  BY  SERIES  FORMULATION 


The  series  type  approach  for  determination  of  vorticity  is  a  convenient 
way  of  satisfying  the  general  upwash  equations. 

V  1.  Series  Forms  for  Trailing  Vorticities  on  the  Wing  and  the  Slat 

We  write  for  the  wing  with  reference  to  Section  III  22splitting  £w  into  two 
contributions  and 

SwC*)  5  SW(©A)  .  Sf  F  S*  ....  (9) 

“here  e*  c-  (.  ^ ) 

If  the  flow  is  attached  at  the  wing  leading  edge  then  the  first 
contribution  is 


Sr  (. ^  -  -  $  t  ^  0A  .  . 

It  is  implied  that  Ti  =0. 

If  the  flow  on  the  wing  is  separated  at  the  leading  edge,  then 
^F  O  and  f,  O 

The  second  contribution  is  denoted  by 

I* 

Sole,)  •  V.  )  ar 

p9l 

The  general  form  for  SW(®A']  is  given  by 


(10) 


(11) 


Sw(.e^  --  -  §,.v.  cef-e^  4  v  y~  <>*(?■ 


(12) 


The  trailing  vorticity  Ss  on  the  starboard  slat  B  comprises  two 
portions  with  reference  to  Section  III  2.2 


^s^s)  -  »  ^&C©b)+  Sb(9b^  ....  (13) 


Mhere  9,  (>/c,*) 

If  the  flow  is  attached  at  the  slat  leading  edge,  then 

Ve»)  m  §3  V.  -W  o.rxd  r3  -  o  ....  (14) 

If  the  flow  is  separated  at  the  slat  leading  edge,  then 
(®»)  -  O 


a.r\d 


r3  ^  o 
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The  second  contribution  is  denoted  by 

s»(A)  -  v-  Y-  b*'  °^Qb) 


The  general  form  for  Se  is  given  by 


.  .  .  .  (15) 


SjCeg.)  ^  V.  gs  fiux  9%  *  V  Y_  H 


.  .  .  (16) 


The  expressions  for  trailing  vorticity  &p  on  the  port  slat  C  follow 
from  geometry  considerations,  i.e. 


Sp  (  ft)  -  -  (^) 

where  Q  is  the  image  of  point  (9-  in  x  z -plane. 

C  6 

<9t  -  — "  (>/c,0 

V  2.  Series  Forms  for  Zero  Total  Circulation  on  the  Slat 
(a)  Attached  flow  about  slat  L.E. 

Equation  (3)  becomes 


.  .  .  .  (17) 


rz  + 


dy 


On  substitution  of  ^  »  -c,x.c«r5S 


we  obtain 


f  +' 

r»-  J  *(.-£)■«<(-£)  -  W  ^ 

J .  I 

n  r 11 


Substituting  Equations  (14)  and  (15) 


r2  +  q  x  Jz  .  V  J  (a^  gp  ^  ©B  +-  V j  hr  •  7  /M*®8)'>w>(9B  o 


^*1  t> 


q  f  V.  3,  £,* -IT  ♦  V.  c.x  .b  £  =  o 
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In  non-dimensional  terms 


r*.  - 

■ — r  +-  J-.  • 

v.  k  *  1 


C|  b  TT  _  r\ 

T  '  '  2  -  ° 


(b)  Separated  flow  about  slat  L.E. 
Equation  (^)  becomes 

C,x 

^  +r3  +  /  Ss 


-c 


(X 


=  O 


which  on  using  similar  analysis  as  in  the  previous  case  yields 

+  r3  +  q  x  .v  b,.  7^  =_  o 


In  non-dimensional  terms 


%  hi  Jk  o 


V.k: 


'  S  C, 

V.K  ^  k 


Equations  (13)  and  (19)  therefore  offer  convenient  relationships  to  remove 
one  of  the  unknowns,  say  bj,  from  the  calculations. 

V  3.  Series  Forms  for  Bound  Vorticities 

The  expressions  for  bound  vorticities  follow  from  application  of  the 
equation  of  conservation  of  vorticity,  i.e. 

dS  .  3*  a 

x —  T  -r —  —  vj  •  ■  * 

<3*  ^ 

Wing 


On  the  wing  we  have 

C*r>  9  ~  \ 

**  9*  Tx  =  T  x* 


i.e. 


and 


nvr\ 


% 


ci6p 


X 

k  x. 


Equation  (20)  becomes 

'  /  st#  -  /*  r©.1  *kc”  *sv  •  • 


We  split  into  three  contributions  as  follows 

rw(6>,i  -  rP  -*■  rs(s,)  ♦  dr 


(18) 


(19) 


(20) 


(21) 


.  .  (22) 


corresponds  to  £f(<9a)  and  arises  only  when  the  flow  is 
attached  at  the  wing  leading  edge  .  V8(fiA)  corresponds  to 

arises  only  when  the  flow  is  separated  at  the  wing  leading  edge. 

On  partial  differentiation  with  respect  to  ,  equation  (22)  yields 


()  <9  A 

Derivation  of  yF  (e^ 


(^Vr  (6fk) 


Using  Equations  (10)  and  (21) 


_  k  Or,  eA  ■  3,  V  £* 


which  gives  by  integration  with  respect  to  within  limits  of  0  and  an 
expression  for  yF  .  * 


YrW  .  -V  9,  v  l  ^ 9 


k  5,  ■ »//’“». 


We  may  also  write  this  in  form 


(-  L ' 


k  3,  V 


kx  / 


which  corresponds  to  linear  attached  flow  solution  of  Jones  slender  wing  approach. 

Derivation  of 

Using  Equations  (15)  and  (21) 


dfcCfc) 


_k  ^  V.  ^ 2F  ^ 

^  **”  \ 
_k.V.  f>.  CLp  ■j'  Co(zM-')e-t  C^Op-O&'j 


which  gives  by  integration  with  respect  to  <9*  within  limits  of  0  and  <9A  an 
expression  for 

.....  .  .  _  J  oU,(2|>-0<9.  ? 


ys(M  .  _k.v. 


L 


Zp-rl 


l 

J 
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Slat  B 

To  obtain  expressions  for  Ys  from  £>s  we  need  to  apply  the  equation  of 
conservation  of  vorticity  (equation  20))  in  terms  of  co-ordinates  x,  yQ,  zQ 
as  in  Figure  20.  This  approach  then  allows  for  inclination  of  the  slat 


%  -  (f*  W  +  oilv  _  c,  <9S  } 


or 


Coo  (9 


8 


% 

C..X 


i.e.  £L 

s  31  = 

Equation  (21)  becomes 


JL 

c-,*1 


Cb  Co,  <7 

and 


~  1 <3) 

Cos  (9 

Cs,®»  35.  • 


3>ics>> 

Te„~ 


A  A 


s 


My  35^(6*)  ✓  . 

*  b  »  (c6Cos«jp  +  ^n^Cp- C,  C*»£g) 


c)  $s(<9b) 
d<9* 


Derivation  of 


Using  Equations  (14)  and  (26) 

_  (Co  to’Cf  4-  Ko^f  -  C>  )  5.  V .  j-  <ySc^t~) 

Z9a 

-  5s  v  t  {  Koou.4?)  oec^)-  «,( 


9  _  /-ICC 


which  gives  on  integration  with  respect  to  $3 

v  *C»  )^(x)  -  .  .  • 

where  F  is  a  constant  of  integration  and  it  is  evaluated  by  using  the  Field 
condition  (equations  (3,  4). 

At  0g  =  0  we  have 

Yc  ( 9s  * 0 )  r  ^  .... 

Derivation  of  Vg  (.^S  ) 


Using  Equations  (15)  and  (26) 
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(28) 


-  19  - 


£>®s 


n 

*  ^  •  C  Cayrf  +  b0ov*  «p  -  c(  ^  ^  •  V '  ^  ®a 

N  V7 

?  V  ^  ^  •  vj^0^  4K0ou;f  )  Gp^  -  CfC*^  c»'V0ffl 

tysl  —i 

N  p 

c  y  (CoU»f  +  K^ui^o.  <y©e 

^:l  - 

-  Si  |"  c«>  (<V'H)®b  +•  '’‘^(ty'l  )(JB 


which  gives  on  integration  with  respect  to  0 

n  r 

V.  ^_*>V  *V^f 


V 


®. 


r  ^Cvo^ .  ^r<v->)g«  \ 
4  {  ~**i  I/*' 


+V.F 


B 


vortex 

"feeding" 

cuts  and 

M 

•I 

0 

It 

drx 
'  Tx 

i.e. 

<5f l 

F*  * 

-*  rB 

(291 


where  F^  is  a  constant  of  integration  and  it  is  determined  from  boundary 
conditions  at  the  edges  of  the  slat.  It  arises  from  q  =  1  term  only  when 

sin(q-l)0R 

- - - reduces  to  0D. 

q-1  0 

We  examine  q  =  1  term 

=  V.  b,  [(Co^cf  +  K0,^)^fi8  -  <k 

If  the  flow  is  separated  at  both  edges  of  the  slat  then  we  allow  for  two 


(30M 


and 


.  -v.  b,  \  it  +  v.  r 


2u 


We  now  have  a  relationship  between 
T3,  f2  and  bi,  i.e. 


d«i 

.  vt, 

_ 

i 

tg 

X 

or 

(dfi  a  £ 1 

‘  ^Fx  *  Jx  ; 

>■  K 

This  relationship  satisfies  the  field  condition  of  zero  total  circulation 
about  the  slat. 


If  the  flow  is  attached  at  the  slat  leading  edge  but  separated  at  its 
trailing  edge,  then  we  allow  only  for  one  vortex  "feeding"  cut  and  write 

Ts,  (<9r  =-rr)  -O 

Equation  (30)  then  yields 


fB.tlfe'T)  ,  O  *  " 


+  V.fB 


Substituting  for  Fg  in  Equation  (30),  we  get 


YBiceB)  r  v.  ^  , 

— 


At  6  =  0,  we  have 
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"^6,  ■-  V.  VCi  TT4_ 


Total  bound  vorticity  at  slat  trailing  edge  9  =  0  is  given  from 

D 

(&sr°)  -  Xi  (  ®6:0)  '+*g(®s:C>) 

*  V  F.  +  v.  Vc,.  1^ 

By  continuity  this  is  fed  into  the  vortex  1*2  using  the  appropriate  sign  at 
the  slat  trailing  edge 

Ts(v») .  -f:  .  -  s 

or  “  &  r  V  F.  +  v.  b,  .  C,.  Tfc 

X  *- 

or  Fv  +■  i  V.  F  V .  C-kZ  .  b,  ■  O 
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Comparing  this  with  Equation  (28),  we  obtain 
F  *  9j  C.  n- 
Y&(6»)  from  Equation  (27)  is 

Y&(Ss)s  3iV-  {  (>«•*  4Ke^+C()tiu^  (V71)!  .  .  .  .  (33) 

V.4  Upwash  Equation  for  Wing 

The  series  forms  for  6W,  6S  and  6p  are  substituted  in  Equation  (1). 

We  distinguish  now  between  the  four  possible  configurations  A-S-A,  S-S-A, 

A-S-S  and  S-S-S. 

(a)  Configuration  A-S-A 

Using  Equations  (1),  (9),  (13) 
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0 

Using  series  froms  for  trailing  vorticity  distributions  from  Equations 
(12),  (16)  and  writing  for  constant  in  the  crossflow  plane.  All 

distances  y  and  z  are  non-dimensionalized  by  dividing  by  kx  • 
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•»,  S'-’iC6*) 


where  the  unknowns 
terms  are  velocity 

(b)  Configuration 

+  (*/. , 
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, l, J *  a  %( ^ 
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(34 


are  9,  ,  ,  bv ,  ?3  and  f*_  .  The  other  multiplying 

influence  coefficients  which  are  defined  in  the  Appendix.  I- 


S-S-A 


+  {  Ss/ 

,^,0)  +  =  o 


Using  series  froms  for  trailing  vorticity  distributions  from  Equations 
(12),  (16),  we  have  for  constant  3^  in  the  crossflow  plane 
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....  (35) 

where  the  unknowns  are  o^> ,  3a  ,  kty ,  H  and  G  .  The  multiplying  terms  are 
the  velocity  influence  coefficients  which  are  defined  in  Appendix  I. 

(c)  Configuration  A-S-S 
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Using  series  forms  for  trailing  vorticity  distributions  as  in  the  previous 
case,  we  have  for  constant  in  the  cross-flow  plane 

M 
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p:\ 

V' 
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....  (36) 

where  the  unknowns  are  H  ^  and  ^  .  The  multiplying  terms 

are  velocity  influence  coefficients,  which  are  defined  in  Appendix  I  . 

(d)  Configuration  S-S-S 
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Using  series  forms  for  trailing  vorticity  distributions  as  in  the  previous 
case,  we  have  for  constant  in  the  cross-flow  plane. 
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where  the  unknowns  are  ^p,  ^  ,  H_  and  fj  .  The  multiplying  terms 

are  velocity  influence  coefficients  which  are  defined  in  Appendix  I. 


(37) 


V  5.  Upwash  Equation  for  Slat 

The  series  forms  for  ,  ^s.  and  are  substituted  in  Equation  (2). 
We  distinguish  as  in  Section  V.4  between  the  four  possible  configurations 
A-S-A ,  S-S-A ,  A-S-S  and  S-S-S. 


(a)  Configuration  A-S-A 

Using  Equations  (2),  (9),  (13) 
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Using  series  forms  for  trailing  vorticity  distributions  from  Equations  (12) 
and  (16)  and  writing  for  constant  in  the  cross-flow  plane 
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The  unknowns  are  3t  ,  ap  ,  5j  ,  ^<y  and  (j  .  The  other  multiplying  terms 
are  velocity  and  the  influence  coefficients  which  are  defined  in  Appendix  I. 

(b)  Configuration  S-S-A 
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Using  series  forms  for  trailing  vorticity  distributions  as  in  the  previous 
case,  we  have  for  constant  in  the  cross-flow  plane 


{  n. C t a.. . )**  -  J 
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^  r" 
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+  vX  f  3°^  (%i- -3^  (*»,%.  >“4  } 

=  —  /OUv,  ....  (39) 

The  unknowns  are  <lp  ,  ky  ,  3j  ,  *7  and  rj_  .  The  other  multiplying  terms 
are  the  velocity  influence  coefficients  defined  in  Appendix  I. 

(c)  Configuration  A-S-S 
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Using  series  forms  for  trailing  vorticity  distributions  as  in  the  previous 
case,  we  have  for  constant  in  the  cross-flow  plane 
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The  unknowns  are  frt  »  *!>  .  ^  ^  and  .  The  other  multiplying  terms 

are  velocity  influence  coefficients  which  are  defined  in  Appendix  I. 

(d)  Configuration  S-S-S 
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Using  series  forms  for  trailing  vorticity  distributions  as  in  the  previous 
case,  we  have  for  constant  in  the  cross-flow  plane 


M 

P~  \ 


*  kw5,(  If,  »*/?««. 

4  lvt«  (*V,  V  .?sJ j 


+  p-  j  3^  ( 

V.k*  ( 


3^  (S/b.Jm)"-1?} 
8°}^  ( K*s 

9^  ( ifce.wW] 


-  -  .  .  .  . 

The  unknowns  are  C*-p,  ,  'i  ,  and  fj  .  The  other  multiplying  terms 

are  velocity  influence  coefficients  which  are  derived  in  Appendix  I. 

(e)  Slats  Only  Configurations  S-A  and  S-S 

The  upwash  equations  can  be  derived  from  configurations  A-S-A  (or 
S-S-A)  and  A-S-S  (or  S-S-S)  by  omitting  the  wing  effects. 


V  6.  Combined  Upwash  Equations 

We  now  incorporate  the  conditions  of  zero  total  circulation  on  the 
slat  (from  Section  V.2)  into  the  upwash  equations  for  the  wing  and  slat 
(from  Sections  V.4  and  V.5). 

Configuration  A-S-A 

(a)  Wing 

Using  Equation  (18)  we  can  eliminate  bj  from  Equation  (34),  i.e. 
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(b)  Slat 


Using  Equation  (18),  bj  can  be  eliminated  Sow  Equation  (38). 
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We  may  write  equations  (42)  and  (43)  in  shorthand  form 

&>‘i'  +  Xtt'p'*r  *  2jA  *  /L 


where  the  unknowns  are 


Ji  ,  «-f  (<’*»•••  M)  >  >  *v  (‘*"i  • '  ^  -  fL 

and  the  coefficients  h  ,  *P  ,  ^5  ,  ,  dt  ,  t-  have  the  obvious  significance 

with  reference  to  equations  (43)  and  (44). 


Configuration  S-S-A 


(a)  Wine 


Using  equation  (18),  bj  can  be  eliminated  from  equation  (35),  i.e. 
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(b)  Slat 


Using  equation  (18),  bj  is  eliminated  from  equation  (39). 
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Equations  (45)  and  (46)  may  be  written  in  shorthand  notation. 
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where  the  unknowns  are 

M),7S,  **(*»*••  *0  ,  ^  '  vTk* 

and  the  coefficients  ,  ?i .  A  -  *>  ’  have  the  obvious 

significance  with  reference  to  equation  (45)  and  (46). 

Configuration  A-S-S 

(a)  Wing 

Using  equation  (19),  bj  can  be  eliminated  from  equation  (36),  i.e. 
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(b)  Slat 


Using  equation  (19)  bj  can  be  eliminated  from  equation  (40),  i.e. 
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Equations  (49)  and  (49)  may  be  written  in  shorthand  notation 
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where  the  unknowns  are  , 

.  and  the  coefficients 
significance  with  reference 
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*  eP  »  have  the  obvious 

to  equations  (48)  and  (49). 


Configuration  S-S-S 


(a)  Wing 


Using  Equation  (19),  b:  can  be  eliminated  from  equation  (37). 
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(b)  Slat 


Using  Equation  (19)  bj  can  be  eliminated  from  (41), 
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Equations  (51)  and  (52)  may  be  written  in  shorthand  notation 
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where  the  unknowns  are  Q-y  (p  =  1  ...  M)>  ^(q  =  2  ...  N)>  -7-  , 

r  p  f  *  V'  »x 

»  “\Tk-x.  an<*  t^le  coe^^cient:s  £/*»  "ff  »  A-i  t  ,  ^3  have  the 
obvious  significance  with  reference  to  equations  (51)  and  (52). 

Slats  Only  Configurations  S-A  and  S-S 

The  combined  upwash  equations  including  the  relationship  of  Equation  (18) 
for  Configuration  S-A  and  Equation  (19)  for  Configuration  S-S  are  derived  by 
omitting  the  wing  effects  in  Equation  (44)  for  Configuration  S-A  and  Equation 
(47)  for  the  Configuration  S-S. 

V  7.  Solution  Procedure 


The  combined  upwash  equations  (44),  (47),  (50)  and  (53)  for  the  various 

types  of  configurations  are  all  sets  of  linear  simultaneous  equations.  The 

procedure  for  solution  by  collocation  at  a  number  of  points  on  the  surfaces  of 

wing  and  slat  is  illustrated  in  Figure  21.  The  various  steps  are  as  follows: 

(a)  Choose  type  of  configuration  (i.e.  edge  conditions). 

(b)  Choose  wing  and  slat  geometry  parameters. 

(c)  Select  (M  +  1)  points  along  wing  semi-span  (equal  intervals  in  0  ). 

(d)  Select  (N  +  1)  points  along  the  slat  span  (equal  intervals  in  0D). 

D 

(e)  Calculate  wing  and  slat  vorticity  induced  velocity  influence  coefficients. 

(f)  Locate  initial  vortex  positions. 

(g)  Calculate  vortices  induced  velocity  influence  coefficients. 

(h)  Solve  the  appropriate  upwash  equation  to  determine  the  unknowns  -  vorticity 
strength  and  vortex  strength. 

(i)  Evaluate  velocities  at  vortex  positions  for  applying  the  equations  for 
zero  force  on  vortex-cut  arrangements. 

(j)  Estimate  predicted  positions  of  vortices. 

(k)  If  predicted  positions  of  vortices  are  not  within  a  pre-specified  small 
tolerance  on  the  initial  assumed  vortex  positions,  then  the  vortices  are 
re-located.  This  re-location  procedure  is  optimised  for  the  best 
direction  of  movement.  Calculation  is  taken  back  to  step  (g). 


(1)  Bound  vorticity,  lift  distribution,  lift  and  drag  are  calculated  (see 
Appendix  II). 


V  8.  Similarity  Parameters 


In  conical  wing  theory  problems,  the  similarity  parameters  are  generally 
taken  as  follows: 


sin  a 


is  the  angle  of  attack  parameter  which  is  the  main  variable  of  the 
problem  and  defines  a  family  of  solutions  related  to  wing  sweepback  and  a. 


1 2  CL 

The  parameters  for  vortex  strength  and  —  for  lift  are  dependent 


on  sin  a/k« 


For  wing-slat  configurations  also,  these  parameters  are  applicable.  The 
angle  of  attack  parameter  that  is  more  appropriate  is  sin  a/k3  referred  to 
slat  leading  edge.  A  family  of  solutions  is  therefore  implied  in  each 
calculation  for  given  values  of  a  and  k3.  Additional  geometry  parameters 
defining  the  gap  between  the  wing  and  slat  may  also  be  identified.  An  under¬ 
standing  of  these  will  lead  to  reducing  the  computations  required. 


V  9.  Computer  Calculations 

Computations  were  performed  on  the  University  of  Bristol  system  4-75 
computer.  The  time  taken  for  each  case  depended  on  the  number  of  collo¬ 
cation  points  on  the  wing-slat  geometry  and  the  accuracy  of  starting 
positions  of  the  vortices.  In  the  order  of  numerical  computation  time 
required  to  attain  reasonable  convergence,  the  various  configurations 
stand  as  follows: - 


(i) 

(ii) 

( iii) 

(iv) 

(v) 

(vi) 


Slats  only  configuration 
Slats  only  configuration 
Wing-slat  configuration 
Wing-slat  configuration 
Wing-slat  configuration 
Wing-slat  configuration 


S  -  A 

S  -  S 

A  -  S  -  A 

S  -  S  -  A 

A  -  S  -  S 

S  -  S  -  S 


It  must  be  mentioned  that  for  configurations  A  -  S  -  S  and  S  -  S  -  S, 
the  convergence  rates  were  rather  slow.  An  improved  iteration  procedure 
may  be  required  based  on  "nested"  iterations  for  each  vortex  in  turn.  A 
different  form  for  series  distributions  of  vorticity  can  also  be  used.  These 
may  have  a  favourable  influence  on  the  convergence  rates  (see  Reference  7). 


40 


VI  RESULTS 


An  understanding  of  the  slats  only  configuration  is  an  essential  prelude 
to  the  wing-slat  configurations. 

VI. 1  SLATS  ONLY  CONFIGURATIONS 


A  series  of  slat  geometries  S1-S3  as  illustrated  in  Figure  22  have  been 
used.  The  slat  leading  edge  has  been  defined  by  y  =  0.25  x  (angle  of  sweep- 
back  75.964  ).  The  trailing  edges  are  given  by  y  =  0.05  x,  0.1  x  and  0.2  x. 

This  forms  a  sequence  with  the  limiting  case  of  the  wing  only  configuration. 

Configuration  S-A  refers  to  the  attached  flow  on  the  slat  leading  edge 
whilst  in  S-S,  flow  separation  at  the  leading  edge  is  included  (see  Figure  16). 

VI  1.1  Configuration  S-A 

The  effects  of  angle  of  attack  ct  on  slat  combination  S1-S2  are  illustrated 
through  Figures  23  and  24.  The  variations  of  load  distributions,  vortex  height, 
vortex  strength  and  lift  coefficient  C^  have  been  considered.  The  effects  of 
increasing  a  are  as  follows 

(i)  The  magnitude  of  slat  wake  vortex  strength  increases  linearly. 

(ii)  The  vortex  moves  inwards  and  outwards  starting  from  a  spanwise  location 
y 2  =  2  (co-cl)  x»  z2  =  0  at  a  =  0. 

(iii)  CLs  increases  linearly. 

The  effect  of  slat  span  on  load  distributions  and  wake  vortex  at  a  =  5° 

O  o  r 

and  10  for  4>  =  0  is  illustrated  in  Figure  25  (a)  and  (b).  This  is  con¬ 
veniently  interpreted  as  the  effect  of  introducing  a  symmetrical  concial  gap 
near  the  centreline  of  a  flat  delta  wing.  The  load  distributions  therefore 
follow  a  plausible  sequence. 

Movement  of  the  trailing  edge  outwards,  i.e.  increasing  the  gap  in  the 
wing  leads  to  the  following  effects  at  fixed  a  (see  Figs  26  and  27). 

(i)  Vortex  strength  magnitude  ^v. (.CptC, reduces  .  However  vortex  strength 
magnitude  ^/^v.  ic,  x  ^  based  on  slat  span  increases  slightly  and 

then  reduces. 

(ii)  Vortex  position  moves  outwards  and  downwards,  the  limiting  values 

being  for  no  gap  and  ilI05  ,  +c 

for  vanishing  slat  span.  '  1  -* 

(iii)  the  total  lift  decreases,  however,  lift  coefficient  based  on  local 
span  increases  slightly  and  then  tends  to  a  limiting  value 
2TT.okw>«C  c©>  A-  (for  an  infinite  sheared  wing  of  sweepback  angleA.  ). 

The  effect  of  slat  inclination  <p  is  shown  for  slat  combination  SI  and  S2 
in  Figures  28  and  29. 

The  effective  angle  of  attack  increases  with  $  increasing  and  aero¬ 
dynamic  effects  may  be  interpreted  on  this  basis,  e.g. 


(i)  The  vortex  strength  magnitude  and  lift  increase  with  $. 
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(il)  The  vortex  spanwise  position  moves  outwards  for  both  increasing  and 

decreasing  <p.  Vortex  height  increases  with  reducing  $  and  conversely. 

(iii)  On  lift  coefficient  basis  (based  on  actual  area  and  not  the  projected 
area),  Cls^qd  improves  with  decreasing  (i.e.  same  sense  as  leading 

edge  droop).  It  must  be  mentioned  that  leading  edge  suction  effect  on 
induced  drag  Cp^  has  not  been  included. 


VI  1.2  Configuration  S-S 

A  few  selected  combinations  S1-S3  from  Figure  22  have  been  used. 

The  effect  of  slat  size  on  load  distributions  and  vortices  is  depicted 
in  Figure  30.  This  is  conveniently  interpreted  as  the  effect  of  introducing 
a  symmetrical  conical  gap  near  the  centre-line  of  the  flat  delta  wing.  It 
is  interesting  to  note  that  peak  loads  are  of  the  same  order,  i.e.  leading 
edge  vortex  lift  dominates  and  this  is  not  substantially  altered  by  the 
vicinity  of  the  slat  trailing  edge. 

The  vortex  strength  and  are  shown  in  Figure  31.  The  magnitude 
of  m.  reduces  with  reduction  in  slat  span  but  tends  to  a  constant  value  for 
vanishing  span.  reaches  a  peak  for  <^02  and  then  decreases 

with  increasing  but  tends  to  constant  value  for  vanishing  span. 

Lift  coefficient  based  on  fCo-Ki)*  decreases  with  increasing 
(Figure  32)  and  tends  to  a  constant  value  for  vanishing  slat  size.  Based 
on  local  span  zc,x  the  curve  as  shown  in  Figure  32  is  obtained,  which 
tends  to  infinity  for  vanishing  slat  size.  This  behaviour  for  vanishing 
slat  size  appears  to  be  reasonable  since  the  usual  incidence  parameter 
sin  “/c,  for  slat  will  also  tend  to  infinity  and  non-linear  lift  is  a 
direct  function  of  this  incidence  parameter. 


The  effect  of  angle  of  attack  a  =  2.5°  to  8°  on  load  distributions  and 
the  vortices  for  slats  geometry  SI  is  shown  in  Figure  33.  The  peak  load 
moves  inwards  with  increasing  o.  ?2  and  T3  both  move  inwards  and  upwards 
with  increasing  o.  Vortex  height  for  all  geometries  increases  with  a  increasing. 

For  a=5°  case,  a  comparison  of  configuration  S-A  and  S-S  is  presented 
in  Figure  34  (i.e.  this  is  the  effect  of  including  T3  at  the  slat  leading 
edge).  We  note  that  with  T3  present,  large  increases  in  lift  distribution 
(non-linear  lift)  and  hence  wake  vortex  strength  are  obtained. 
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It  is  of  interest  also  to  compare  the  vortex  and  lift  characteristics 
for  the  two  configurations  S-A  and  S-S  throughout  the  incidence  range 
(Figures  35  and  36).  We  note  that  by  including  Ts,  large  benefits  in 
non-linear  lift  are  obtained.  13  approaches  the  strength  given  by  Brown 
and  Michael  (for  the  wing  with  the  same  leading  edge).  This  is  particularly 
interesting  and  shows  that  the  effect  of  slat  trailing  edge  on  T 3  and  hence 
on  non-linear  lift  is  small. 

The  magnitude  of  for  the  configuration  S-S  is  higher  and  increases 
non-linearly  with  a.  This  is  in  line  with  total  zero  circulation  condition 
about  the  slat. 

An  idea  of  the  /Cp  characteristics  can  be  gained  from  Figure  37. 

As  the  slat  span  reduces  (i.e.  it  approaches  the  case  of  infinite  sheared 
wing),  its  performance  increases,  as  might  be  expected.  As  the  slat  gap 
vanishes,  the  limit  is  the  Brown  and  Michael  curve. 

The  comparison  of  vortex  locations  of  and  T3  presented  in  Figure  38 
refers  to  the  case  of  twin  slats  (combination  S3)  compared  with  the 
case  of  one  slat  only  (i.e.  yawed  delta  wing  results  from  Cohen^3-3)).  It 
is  noted  that  although  ?3  locations  are  comparable,  ?2  location  for  the  twin 
slats  combination  is  higher  because  of  the  mutual  influence  of  one  slat  on 
the  other. 
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VI  2.  WING-SLAT  CONFIGURATIONS 

A  series  of  wing-slat  geometries  WS1-WS20,  as  shown  in  Table  I  ,  have 
been  used  for  analysis  and  comparison.  In  most  cases  the  slat  leading  edge 
(at  4>  =  0  )  is  at  y  =  k3X  =  0.25x  and  the  maximum  angle  of  attack  a  is  of 
the  order  of  15°.  This  gives  the  maximum  value  for  angle  of  attack  similarity 
parameter  (sin  a/k3>  as  1.034. 

There  are  four  configurations  of  Figure  15  to  be  considered  (A-S-A, 

A-S-S,  S-S-A  and  S-S-S)  all  with  different  combinations  of  edge  conditions. 

VI  2.1  Configuration  A-S-A 

The  conical  streamline  pattern  sketched  in  Figure  39  results  from  the 
effect  of  introducing  a  conical  gap  in  a  flat  delta  wing  away  from  the  wing 
centre-line  (constituting  thus  a  wing-slat  configuration).  Load  distributions 
for  a  typical  case  at  a  =  5  are  shown  in  Figure  40.  These  have  been  compared 
with  the  cases  of  slats  only  (Configuration  S-A)  and  wing  only  (no  gap  - 
Jones'  linear  theory  result).  It  is  evident  that  wing  and  slat  in  combination 
produce  more  lift  than  the  sum  of  their  individual  lift  values,  but  it  is  less 
than  the  value  for  the  wing  only  as  might  be  anticipated.  The  appearance  of 
non-linear  lift  on  the  wing  under  the  slat  trailing  edge  wake  vortex  is  also 
noted. 

The  effect  of  angle  of  attack  a.  =  2.5°  toll?  on  load  distributions  and 
vortex  position  is  shown  in  Figure  41.  The  slat  wake  vortex  moves  inwards 
and  upwards  with  increasing  a  as  its  magnitude  increases  (Figure  42).  It  is 
interesting  to  compare  this  strength  with  the  equivalent  value  from  slats 
only  conf iguration  S-A  and  it  is  noted  that  the  presence  of  the  wing  leads 
to  an  increase  in  the  vortex  strength  magnitude. 

Lift  coefficient,  angle  of  attack  relationships  are  shown  in  Figure  43. 

It  is  convenient  to  distinguish  between  the  lift  coefficients  based  on  the 
total  projected  area  (between  the  slat  leading  edges  at  <p  =  0  )  and  the  actual 
exposed  surface  area  (i.e.  allowing  for  the  slot).  It  is  noted  that  the  wing 
slat  case  offers  lift  of  the  same  order  as  the  wing  with  the  same  exposed 
surface  area. 


The  effect  of  slat  inclination  =  0  to  -  20  (in  the  same  sense  as  wing 
leading  edge  droop)  is  shown  in  Figure  44.  It  is  noted  that  the  'ip  between 
the  wing  leading  edge  and  the  slat  trailing  edge  increases  as  <t>  de  -es. 

The  lift  distribution  on  the  slat  and  hence  the  magnitude  of  the  s_  wake 
vortex  decreases  with  decrease  in  Similarly  the  effect  of  slat  wane  voi  i  x 
on  the  wing  load  distribution  also  decreases  as  d>  decreases. 

The  effect  of  gap  size  for  various  <p  has  also  been  examined  (Figure  45a, 
b,  c).  Widening  of  the  gap  leads  to  loss  of  lift  distribution  on  the  slat  as 
may  be  expected  and  this  leads  to  a  reduction  in  non-linear  peak  loading 
induced  on  the  wing.  For  the  <f>  =  0  case,  the  limiting  lift  distribution  for 
vanishing  gap  (overall  wing)  case  is  also  illustrated  as  a  reference  curve. 

The  effect  of  the  slat  height  is  depicted  in  Figure  46.  We  note  that  as  the 
slat  distance  reduces,  the  slat  loading  and  its  effect  on  the  wing  increases, 
i.e.  the  loading  near  the  wing  leading  edge  decreases  and  the  peak  load  on  the 
wing  under  the  slat  wake  vortex  increases.  For  negative  slat  height  the  slat 
wake  vortex  lies  nearer  to  the  wing  leading  edge. 
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It  is  of  interest  next  to  consider  the  efficiency  of  wing-slat 
conf igurations . 

Figure  47  shows  the  possible  design  approach  for  a  high  speed  wing. 

It  might,  for  example,  feature: 

(1)  a  hinged  flap; 

(2)  leading  edge  slat  which  can  have  a  variable  gap  and 

deflection ; 

(3)  leading  edge  deflection  and  extension  with  major  geometry 

changes. 

(1)  and  (2)  above  are  comparable  since  there  is  no  change  in  plat¬ 
form  area,  but  (3)  really  means  an  increased  wing  area. 

Results  from  a  selection  of  geometries  (k  =  0.18,  cQ=  0.2  and  Cj  =  - 

have  been  compared  on  Cb/Cr^  and  Cg  basis  with  a  flat  wing  of  the  same  are  i 
in  Figs.  4a  and  49.  The  leading  edge  suction  thrust  term  has  been  discoun¬ 
ted  in  Fig. 48  but  retained  in  Fig. 49.  It  appears  that  for  the  configura¬ 
tion  A-S-A,  the  benefits  of  the  slat  are  dependent  on  the  realisation 
of  leading  edge  suction  thrust.  If  the  suction  thrust  is  fully  realised 
at  the  leading  edges  then  the  wing  only  case  is  superior  throughout  the 
Cp  range.  In  practice,  however,  it  is  difficult  to  achieve  the  idealise: 
case  of  190%  suction  and  a  20°-  -  90'  figure  may  be  more  realisti  . 
Furthermore,  flee  separations  will  occur  and  alter  the  configuration 
assumptions  at  the  leading  edges.  It  is  therefore  of  interest  to  loo) 
at  the  basL  case  of  Fig. 48  without  leading  edge  suction. 

(i)  For  =  0°,  improvements  in  (pressure)  b/..’-.  are  possible 
for  C.  >  .13,  e.g.  If-'  at  C,  =0.3. 

J_  i  1. 

(ii)  With  regard  to  slat  height  at  <}.  =  0°,  there  appears  to  bo 
an  optimum  position  for  the  best  rL/Cpp. 

(iii)  The  '“L/C;  ,  values  are  appreciably  enhanced  (about  If.)  by 

combining  height  hQ  and  inclination  <t>  of  the  slat  (in  the 
sense  of  leading  edge  droop). 

(iv)  There  is  a  fur'ther  8%  gain  in  Cg  if  the  curves  for  wing-slat 
are  based  on  the  total  exposed  area  of  the  wing  and  slat 
(i.e.  discounting  the  gap). 

(v)  More  optimum  configurations  may  be  derived  by  allowing  slat 
span  as  an  additional  variable. 

It  may  be  interpreted  that  gains  in  LL/Cf).  are  due  to  higher 
operating  efficiency  of  the  slat  which  may  operate  at  a  higher  value  of 
incidence  parameter  aB/eg  than  the  comparable  figure  for  the  wing,  i.e.  a/f  , 

wnere  e  is  semi -apex  angle  of  the  wing  =  tan  ^  k 

and  eg  is  semi-apex  angle  of  the  slat. 

In  the  configuration  A-S-A,  the  wing  basic  lift-incidence  realtionship 
is  given  by 

“  2m  k  sin  a  (based  on  k:<) 

whilst  the  slat  lift-incidence  relationship  is  more  akin  to  an  infinite 
apect  ratio  yawed  wing 

CLs  =  2itcosA  .  ag  (based  on  2  c_l  x )  ■ 

For  certain  combinations  of  a,  c,  ag  and  eg  therefore  the  wing  and  slat 
provide  improvements  in  efficiency  compared  with  the  wing  alone. 
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VI  2.2  Configuration  A-S-S 

The  conical  streamline  pattern  for  this  class  of  flow  is  sketched 
in  Figure  50.  Compared  with  configuration  A-S-A,  the  slat  leading  edge 
and  trailing  edge  vortices  are  expected  to  be  stronger  and  this  will 
increase  the  downwash  effect  at  the  wing  leading  edge. 

The  effect  of  angle  of  attack  a  =  2.5  to  6.5°  on  load  distributions 
and  vortex  positions  is  illustrated  in  Figure  51.  Both  slat  vortices 
move  inwards  and  upwards  with  increasing  a.  The  peak  loads  induced  on  the 
wing  and  the  slat  'under  *he  two  vortices  both  increase  as  a  increases. 

Load  distribution  for  a  typical  case  at  a  =  5°  have  been  compared 
with  cases  of  slats  only  (Configuration  S-S)  and  wing  only  (no  gap  -  Brown 
and  Michael  theory  result)  in  Figure  52.  It  is  evident  that  wing  and 
slat  in  combination  produce  more  lift  than  the  sum  of  their  individual 
lift  values,  but  it  is  less  than  the  value  for  the  wing  only  as  might  te 
anticipated.  The  load  near  the  wing  leading  edge  is  reduced  under  the 
strong  downwash  effect  of  the  slat,  but  the  central  portion  of  the  wing 
carries  increased  loading. 

Load  distributions  and  vortex  positions  for  a  =2.5  and  5°  have  been  .-ompa 
with  those  of  configuration  A-S-A  in  Figure  53.  This  demonstrates  the 
effect  of  including  the  vortex  T3  at  the  slat  leading  edge.  The  slat  in 
the  configuration  A-S-S  carries  increased  load  at  the  expense  of  reduced 
load  on  the  wing  near  the  leading  edge.  The  peak  load  near  the  centre  of 
the  wing  for  A-S-S  under  the  stronger  influence  of  Fp  is  also  higher. 

The  vortex  strengths  Tj  and  Fj  for  Configuration  A-S-S  have  been 
compared  with  the  case  of  slats  only  Configuration  S-S  in  Figure  54.  The 
Brown  and  Michael  curve  for  13  has  also  been  shown,  as  well  as  the  curve 
for  ?2  from  Configuration  A-S-A.  It  is  noted  that  at  a  given  a,  the  value 
of  T3  required  for  Configuration  A-S-S  lies  between  the  narrow  band  for  the 
slats  only  Configuration  S-S  and  Brown  and  Michael,  indicating  that  the 
slat  trailing  edge  has  a  small  effect  on  the  slat  leading  edge  vortex. 

The  magnitude  of  for  Configuration  A-S-S  is  greater  than  that  for 
either  Configurations  A-S-A  or  S-S. 

Lift  coefficient,  angle  of  attack  relationships  are  shown  in  Fig. 55. 

It  is  noted  that  the  wing-slat  configuration  A-S-S  offers  an  appreciable 
increase  in  lift  compared  with  the  wing  of  same  area  according  to 
the  Brown  and  Michael  theory.  In  fact  it  approaches  the  Brown  and  Michael 
curve  discounting  the  gap.  Linear  theory  curves  are  shown  as  a  reference. 

The  effect  of  slat  inclination  $  =  0°  to  -10°  (in  the  same  sense  as 
wing  leading  edge  droop)  is  shown  in  Fig.  56.  With  decrease  in  <t>,  the 
lift  distribution  on  the  slat  and  hence  magnitudes  of  both  Y^  and  1 3 
decrease.  The  effect  of  slat  wake  vortex  on  the  wing  load  distribution 
also  decreases. 

C  C 

Results  from  a  few  calculations  have  been  compared  on  L/C^.  and  L 
basis  with  a  flat  wing  of  the  same  area  in  Fig. 57.  Leading  edge  suction 
thrust  at  the  wing  leading  edge  has  been  included.  It  is  noted  that  the 
wing  slat  configuration  offers  15%  -  20%  gain  in  ^L/Cq^  at  ^L  about  0.3. 
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VI  2.3  Configuration  S-S-A 

The  conical  streamline  pattern  for  this  class  of  flow  is  shown  in 
Figure  58.  Although  the  wing  leading  edge  vortex  is  shown  above  the 
wing  surface,  it  is  possible  for  the  vortex  to  appear  under  the  wing 
surface  because  of  the  strong  downwash  effect  of  the  slat. 

The  effect  of  angle  of  attack  a  =  2.5°  to  15°  on  load  distributions 
and  vortex  positions  is  illustrated  in  Figure  59.  Both  vortices  Fj  and  r2 
move  inwards  and  upwards  with  increasing  a.  The  peak  loads  induced  on 
the  wing  under  the  two  vortices  both  increase  as  a  increases.  It  is  of 
interest  to  compare  the  two  vortex  strengths  with  corresponding  values 
from  slats  only  configuration  S-A  and  Brown  and  Michael  wing  only  case 
(Figure  60).  The  presence  of  wing  leads  to  an  increase  in  magnitude  of  the 
slat  wake  vortex  r2  as  also  noted  in  configuration  A-S-A.  The  slats 
cause  the  wing  vortex  T j  to  decrease  in  strength. 

Load  distributions  and  vortex  positions  for  a  =  2.5°  and  5°  have  been 
compared  with  those  of  configurations  A-S-A  and  A-S-S  in  Fig.  61(a)  and  (b). 
Configurations  S-S-A  and  A-S-A  are  more  readily  comparable  because  slat 
loadings  are  similar  and  the  difference  can  be  interpreted  as  the  effect  due 
to  the  inclusion  of  wing  leading  edge  vortex  Tj.  The  presence  of  Tj  implies 
therefore  an  increase  in  loading  near  the  wing  leading  edge.  In  configura¬ 
tion  A-S-S  the  slat  loadings  are  higher  and  therefore  the  wing  leading  edge 
loads  are  much  reduced. 

Lift  coefficient,  angle  of  attack  relationships  are  shown  in  Fig.  62. 

It  is  noted  that  the  wing-slat  configuration  S-S-A  offers  an  appreciable 
increase  in  lift  compared  with  linear  theory  based  on  both  types  of 
reference  areas  (i.e.  with  and  without  the  inclusion  of  the  slat  gap). 

Results  from  a  few  calculations  have  been  compared  on  ^L/Cp^  and  Cp 
basis  with  a  flat  wing  of  the  same  area  in  Figure  63.  Curves  with  and 
without  leading  edge  suction  thrust  have  been  shown;  Brown  and  Michael 
curve  is  also  shown  for  reference.  It  is  noted  that 

(i)  with  realisation  of  100%  leading  edge  suction  the  configuration  S-S-A 
is  not  superior  to  the  wing  only  case  with  attached  flow,  but  it  is 
superior  to  the  Brown  and  Michael  curve  at  Cp  >  0.3. 

(ii)  If  the  leading  edge  suction  thrust  is  not  included,  then  there  is  an 
appreciable  gain  in  CL/Cpp  at  higher  values  of  Cp  (e.g.  40%  at 
Cp  >  0. 5) . 
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IV  2.4  Configuration  S-S-S 

The  conical  streamline  pattern  for  this  class  of  flow  is  shown  in 
Figure  64.  Although  the  wing  leading  edge  vortex  is  shown  above  the 
wing  surface,  it  is  possible  for  it  to  appear  under  the  wing  surface  due 
to  the  strong  downwash  effect  of  the  slat.  Because  of  the  highly  non¬ 
linear  nature  of  the  flow,  convergence  of  the  solution,  i.e.  the  three 
vortices,  was  found  to  be  extremely  slow  and  only  a  few  examples  have 
been  attempted. 

Figure  65  shows  the  load  distribution  and  vortex  positions  at  a  =  2.5°. 
This  has  also  been  compared  with  configuration  A-S-S  to  show  the  effect  of 
including  wing  leading  edge  vortex  T j .  The  presence  of  T j  appears  to  imply 
a  reduction  in  loading  near  the  wing  leading  edge  and  an  inward  and  upward 
movement  of  T3  with  an  appropriate  movement  of  suction  peak  loading  on  the 
slat.  It  must  be  emphasised  that  since  these  effects  are  highly  non-linear 
and  interdependent,  behaviour  cannot  be  easily  extrapolated  to  other 
geometries. 

The  effect  of  slat  height  is  depicted  in  Figure  66  for  a  =  5°.  Results 
for  smaller  slat  heights  were  found  very  difficult  to  obtain  and  therefore 
only  two  heights  have  been  considered.  It  is  noted  that  as  the  slat  height 
reduces,  the  slat  loading  and  its  effect  on  the  wing  increases,  i.e.  the 
loading  near  the  wing  leading  edge  decreases  and  the  peak  load  on  the  wing 
under  the  slat  wake  vortex  increases. 

It  is  interesting  to  look  at  the  significance  of  leading  edge  conditions 
for  a  given  slat  height.  Figure  67  shows  a  comparison  of  configuration  A-S-A 
with  S-S-S  and  we  note  the  large  changes  in  lift  distribution  near  the  wing 
leading  edge  (due  to  Tj)  and  the  slat  leading  edge  (due  to  Tj).  The  slat 
trailing  edge  vortex  F2  is  increased  in  strength  for  configuration  S-S-S 
and  therefore  causes  a  larger  peak  loading  on  the  wing. 

Results  from  a  few  calculations  have  been  compared  in  Figure  68  on  j 
CL/Cp^  and  CL  basis  with  a  flat  wing  according  to  Brown  and  Michael  theory 
for  both  types  of  reference  areas  (i.e.  with  and  without  the  inclusion  of 
the  slat  gap).  These  results  indicate  a  very  appreciable  30-40%  gain  in 
CL/Cp£  at  CL  values  above  0.2. 
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VII  CONCLUSIONS 


The  conical  type  approach  gives  an  insight  into  flows  around  swept- 
back  slender  wings  with  leading  edge  devices-slats.  The  prevailing 
conditions  at  the  edges  (flow  attached  or  separated)  of  the  configurations 
have  a  marked  effect  on  the  performance.  Separated  flow  conditions 
generally  lead  to  appreciable  benefits  in  lift-drag  ratio,  e.g.  30-40% 
at  higher  lift  coefficients. 

The  method  of  this  report  allows  a  large  number  of  possibilities 
for  variation  of  geometry,  and  further  optimisation  studies  can  be  under¬ 
taken.  From  some  of  the  cases  considered,  it  can  be  inferred  that 
theoretical  gains  in  performance  are  of  the  same  order  as  those  measured 
in  experiments. 

The  leading  edge  devices-slats  can,  therefore,  be  useful  not  only 
in  the  obvious  case  where  high  lift  is  required,  e.g.  near  an  airfield, 
but  also  for  transonic  manoeuvrability  and  increased  performance  and 
controllability  at  both  low  and  high  speeds. 
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VIII  FUTURE  WORK  AND  RECOMMENDATIONS 


There  are  several  areas  of  work,  both  theoretical  and  experimental, 
that  emerge  from  the  present  report  and  it  appears  that  there  is  some  way 
to  go  before  configurations  such  as  those  of  Figures  6  and  69,  Reference 
14,  can  be  tackled  successfully. 

(a)  Theoretical 


The  slender  wing  methods  may  be  developed  as  follows: - 

(i)  Incorporation  of  wing  camber.  An  approach  using  conformal  trans¬ 
formations  is  illustrated  in  Figure  70.  The  method  mentioned  by 
Nangia^-5)  may  ^  used  to  develop  realistic  leading  edge  camber. 

In  addition  to  vorticity  distributions,  source  distributions  will 
also  be  required.  Calculations  are  performed  in  the  simpler 
transformed  plane. 

(ii)  Incorporation  of  wing  and  slat  thickness  and  camber.  Multiple  aero¬ 
foil  type  conformal  transformation  methods  (Ref.  16,  17)  can  be 
used  to  simplify  the  thickness  problems  (Figure  71). 

(iii)  Vortex  sheet  representation  of  separations.  This  is  an  important 
aspect  and  has  a  significant  bearing  on  the  flow  near  the  leading 
edges.  The  work  of  Pullind°>il)  and  Jones^2)  £s  Df  interest.  The 
first  case  to  consider  in  the  inclusion  of  vortex  sheet  at  the  slat 
trailing  edge. 

(iv)  Multiple  slats  are  a  fairly  common  feature  of  conventional  wings  of 
large  aspect  ratio.  They  offer  large  gains  in  lift.  A  theoretical 
study  within  the  framework  of  slender  wing  theory  can  be  formulated 
to  assess  their  potential  (Figure  72).  Conformal  transformation 
method  aids  in  simplifying  the  geometry. 

(18  ) 

(v)  Incorporation  of  planform  effects.  Step-by-step  methods,  e.g.  Smith 
and  Clarke^^)  for  wings  only  may  be  extended  (Figure  73).  This  will 
be  the  first  approximation  to  the  treatment  of  planform  exactly  by 
inclusion  of  chordwise  terms,  e.g.  as  in  methods  of  Nangia^)  and 
Tinoco  and  Yoshihara^4 ) . 

(vi)  Assessment  of  viscid  phenomenon  and  interactions  -  this  is  obviously 
a  difficult  area,  but  some  work  with  conical  type  boundary  layers  is 
in  order. 

The  theoretical  programme  should  provide  a  good  understanding  of  the 
flow  in  realistic  wing-slat  configurations  and  also  aid  in  definition  of 
limits  of  parametric  variation  in  experimental  studies. 

(b)  Experimental 

In  view  of  the  fact  that  there  is  really  a  scarcity  of  data  on  slender 
wings  with  leading  edge  devices,  a  number  of  models  can  be  devised  with 
varying  degrees  of  sophistication  starting  from  conical  type  models. 

(i)  Conical  type  Models  (Figure  74)  enable  an  understanding  of  the  physical 
flow  features  with  relatively  low  cost.  Force  and  pressure  plotting 
models  are  of  interest  here.  Systematic  flow  visualisation  studies 
are  also  required. 
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(ii)  Non-conical  Models  (Figure  75).  A  simple  geometry  that  needs  to 
be  tested  is  the  effect  of  non-conical  slats  on  a  delta  wing. 

The  relative  merits  of  various  types  of  slats  can  then  be  assessed. 

The  effect  of  part  span  slats  on  a  general  planform  is  the 
next  study  required.  It  may  be  argued  that  high  lift  device  is 
best  employed  ahead  of  the  wing  aerodynamic  centre  along  with  a 
flap  type  trimming  surface  near  the  wing  trailing  edge. 

(iii)  Multiple  slats.  An  experimental  study  on  slender  wings  is  required 
to  assess  their  application  potential  and  also  for  evaluation  of 
possible  theoretical  results. 

(iv)  Viscid  and  Interaction  Effects.  Detailed  flow  survey  type  studies 
are  required  to  give  an  idea  of  the  vortex  wake  and  boundary  layer 
interactions . 
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APPENDIX  I 


INDUCED  VELOCITIES  AND  INFLUENCE  COEFFICIENTS 
Wing  Vorticity  5W 

The  velocity  components  V"  and  induced  at  a  general  point  (  ) 

in  the^  cross-flow  plane  are 
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By  substituting  for  the  components  of  Sy, ,  the  velocity 
coefficients  can  be  written  as  follows: 
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where 
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on  the  surface  of  wing. 
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The  expressions  for  9  and  can  also  be  derived  from  potential 

theory  using  complex  variables.  These  are  exact. 
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Starboard  Slat  Vorticity  g4 


The  velocity  components  and  Wfc  induced  at  general  point  ( ^  i 
in  the  cross  flow  plane  (in  slat  axes  system  Jj  ,  *  ).  *ft 
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By  substituting  for  the  components  of  $s  ,  the  velocity  influence  coef¬ 
ficients  can  be  written  as  follows: 
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Vortex  pairs,  ±  r i  ±  F2  ±r3 


The  velocity  components  FT"  and  W"  induced  due  to  a  pair  of  vortices 
±  rt  at  (  )  (i  =  1,  2  or  2),  at  a  general  point  ( 2a  ,  ) 

in  the  cross  flow  plane  are 
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APPENDIX  II 


FORCES  AND  FORCE  DISTRIBUTIONS 

Forces  on  the  Wing 

The  normal  force  coefficient  distribution  on  the  wing  is  given  by 

\l  V  V  J 

where  is  the  total  spanwise  velocity  induced  at  the  point  y .  . 

The  total  normal  force  coefficient  on  the  wing  is  obtained  by 

integration  of  Cj^r  along  the  wing  span.  Lift  coefficient  -  Cp  and  Drag 
coefficient  -  Cp  (not  including  the  profile  drag)  follow  from  W 

ClW  Ctr>< 

° **  ‘X.  +  Q>Wt 

where  C-,  is  the  drag  due  to  leading  edge  suction  on  the  wing  (if  presen 
i.e.  the  «  onf ignrat ions  A-S-A  and  A-S-s) .  Cpw  may  he  calculated  ac:ordi 
to  methods  inti,  iteu  in  Ret.  jo.  It  depends  on  the  value  of  coefficient 

OWv  o<  Tt  k  .  J  I  _  \c.2  .  i 

where  h  rot  cr  •  to  refere;.-.*:  plunlorm  with  semi -span  k  .x 
i e i  ret 

Forces  on  the  Slat 

The  normal  force  coefficient  distribution  on  each  slat  is  given  by 

where  ^b^b)  is  the  total  velocity  induced  in  the  ^  direction  it  the 
point  . 

The  LOtdl  normal  force  coefficient  on  the  slat  Cj^  is  obtained  fy 
integration  of  C';jL  along  the  slat  span. 

Lift  coefficient  (Cp  )  and  drag  coefficient  ((.'p__)  i allow  from 

CL.-  =  cN.3  (:os  as 

(:DS  ~  rHs  sin  as  +  ('r>s. 

The  profile  drag  has  not  been  included. 
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where  Cps  is  the  drag  due  to  leading  edge  suction  on  the  slat  (if  present 
i.e.  in  t^e  configurations  A-S-A  andJ-S-A).  Cps  may  be  derived  on  the 
basis  of  Ref.  14  and  its  value  depends  on  the  coefficient  g3. 

For  one  slat  ,  and  based  on  non-dimensionalising  semi -span  kg  x, 
U'ig.20)  Cq  referred  to  reference  planform  of  semi-span  k^^.x  is  given  by 

CD  =  -  3j.  TT  .  X .  keJT-  1?  .  K- 

Ss  i  S 

vfUere  kg  =  Co  Ctncf  ■+  k0/>w<p  +  ^ 


Total  Forces  on  Wing-Slat  Configuration 

Total  lift  (CL)  and  Induced  Drag  Coefficient  (CD)  are  given  by  adding 
the  wing  and  slat  contributions. 

*“b  ^Lw  +  Lr> 


'ino  reference 


lengths  should  be  consistent. 
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